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Abstract 

Two path integral representations for the T matrix in nonrelativistic potential scattering are de- 
rived and proved to produce the complete Born series when expanded to all orders. They are 
obtained with the help of "phantom" degrees of freedom which take away explicit phases that 
diverge for asymptotic times. In addition, energy conservation is enforced by imposing a Faddeev- 
Popov-like constraint in the velocity path integral. These expressions may be useful for attempts to 
evaluate the path integral in real time and for alternative multiple scattering expansions. Standard 
and novel eikonal-type high-energy approximations and systematic expansions immediately follow. 



I. Introduction 



Nonrelativistic quantum mechanical scattering in a local potential is usually described in the frame- 
work of time-dependent or time-independent solutions of the Schrodinger equation (see, for example, 
Ref. [1]). Path-integral methods in quantum mechanics, on the other hand, are mostly applied to 
the discrete spectrum, e.g., for harmonic [2] or anharmonic oscillators, in particular for evaluating 
the energy splitting in the double- well potential [3]. In contrast, the transition matrix for the con- 
tinuous spectrum is rarely represented as path integral. Even if available, many representations turn 
out to be rather formal, e.g., requiring infinitely many differentiations [4, 5] or infinite time limits 
to be performed [6, 7]. This is not only impractible but also unfortunate since a convenient path 
integral representation may lead to new approximations and may be extended readily to the many- 
body problem or Quantum Field Theory. Also the long-standing problem how to evaluate real-time 
path integrals by stochastic methods needs a suitable path integral representation as starting method. 
There has been significant progress in dealing with real-time path integrals for dissipative systems [8] 
and with coherent-state path integrals for autocorrelation functions [9] but in closed systems and 
infinite scattering times only zero-energy scattering seems to be tractable by Euclidean Monte Carlo 
methods [10] at present (for other attempts, see Refs. [11, 12]. 

Medium- and high-energy many-body scattering has to rely on multiple scattering expansions apart 
from the very restricted few-body cases where exact quantum mechanical calculations are possible 
[13]. One of the most simple and versatile multiple scattering versions - Glauber' s approach - is 
based on the time-honored eikonal approximation where the particle is assumed to travel along a 
straight-line trajectory. This restricts the application usually to high-energy, forward scattering. For 
potential scattering systematic improvements to this approximation have been worked out long ago 
[14, 15] but even at high energies the convergence of these expansions is unsatisfactory for some 
classes of potentials. There are numerous other studies which try to extend the range of validity 
of the eikonal approximation (see, e.g. Ref. [16]). Clearly a path integral representation for the T 
matrix which naturally gives rise to these high-energy approximations would be useful both for the 
analytical and (perhaps) for the numerical problems mentioned above. As one of the merits of a 
path integral approach is its direct extension to field theory one may also expect applications in the 
relativistic domain where the usual procedure for an eikonal approximation consists of simplifying 
individual diagrams and resumming them [17]. 

In the present work we will derive a path-integral representation of the T matrix in potential 
scattering which is similar to the one given by Campbell et al. [18] many years ago (see also Refs. 
[19, 20]). However, ours is not a phase-space path integral as developed there but a particular path 
integral over velocities which is a significant reduction in complexity. The most obvious application is 
at high energy where a new sequence of high-energy approximations immediately follows. However, 
the main aim of the present work is not to give another high-energy approximation but to demonstrate 
that path integral methods lead to new, conceptually (albeit not technically) simple results which may 
be extended to the many-body case. 

Preliminary results have already been presented elsewhere [21] and a (slightly different) account 
is included in a textbook on path integrals [22]. These previous attempts suffered from ambiguities 
in the limit of large scattering times where energy conservation and the elimination of "dangerous" 
phases from the S matrix have to be achieved. In particular, the order in which these procedures 
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were taken seemed to give many, at first sight equivalent formulations which however did not stand up 
to further scrutinity. In the present, detailed account we first eliminate these phases by introducing 
"phantom" degrees of freedom (dynamical variables with the wrong-sign kinetic term) and then isolate 
the variables whose large-time behavior gives rise to energy conservation by a suitable insertion of 
unity into the path integral. This is the classical Faddeev-Popov trick which first was utilized by 
Campbell et al. for path-integral descriptions of potential scattering. The resulting path-integral 
representation of the T matrix is shown to be valid by explicitly working out the Born series to 
arbitrary order. We give two versions of this path-integral representation corresponding to different 
reference paths (straight-line or eikonal and "ray") about which the quantum fluctuations have to 
evaluated. 

This paper is organized as follows. In Sec. II we introduce velocity path integrals which are par- 
ticularly suited for our purposes. Sections III and IV describe how "phantom" degrees of freedom 
naturally arise and the implementation of our particular Faddeev-Popov constraint. The "ray" rep- 
resentation is developed in Sec. V and systematic high-energy expansions are worked out in Sec. VI. 
These are tested numerically for scattering from a Gaussian potential in Sec. VII followed by our 
conclusions and outlook. More technical details can be found in three appendices. 



II. Velocity path integrals for the S matrix 

Consider nonrelativistic scattering in a local potential V(r) which vanishes asymptotically so that the 
corresponding Hamiltonian (also) has a continous spectrum. The initial momentum of the particle 
with mass m is kj (h = 1) and the final momentum kj . Our scattering states are normalized 
according to {4>f\4>i) = (2-7r) 3 (5( 3 )(kf — kj) . Time-dependent scattering is formulated in the interaction 
picture [1] in which the free propagation has been removed. The S matrix is then just the matrix 
element of the time-evolution operator in the interaction picture: 



Ui(t b ,t a ) 



JH t b 



U(t b ,t a )e 



-iHo t a 



U(t b ,t a ) = exp -iH{t b -t a 



(1) 



taken between scattering states and evaluated at asymptotic times 



km U f Uj(T,-T) 

1 — »oo \ 



(Ei+E f )T 

(2n) 3 S^ (kj - k f ) - 2iri5 {Ei - E f ) . 



lim e^^f' 1 {(f> f U(T,-T) 



(2) 



The second line defines the T matrix after the energy conserving 5 function has been factored out. 
Then Ei = k|/(2m) = kj/(2m) = Ef = k 2 /(2m) = E is the common scattering energy. 

To find a path integral representation of the T matrix we start from the standard path integral 
expression for the matrix element of the time-evolution operator U (x b ,t b ;x a ,t a ) = 
< x b \ exp[—iH(t b — i a )]|x a > [2] in which one integrates functionally over all paths starting at x a at 
time t a and ending at x& at time t b . As usual this is realized by dividing the time difference into iV 
intervals e = (t b — t a )/N and integrating over all intermediate points x^ , k = 1, . . . N — 1 with the 
exponential of i times the classical action as weight. 
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For our purposes it is, however, more convenient to integrate functionally over velocities [20, 23] 
which is achieved by multiplying the time sliced path integral for U with the following factor 



N 



i = n / ^ 6 

k=l J 



(3) 



The Xfc integrations can then be performed which gives Xj = xo + e^ = i Vj, or in the continous 
notation the trajectory x(t) = x a + J * (it' v(i') . However, one 5 function remains and we obtain 



U(x b ,t b ;x a ,t a ) 



lim 



, m ,3N I N \ 



•exp<; «e^[yv|-y ( Xj =x a + e^ Vi 



3 

1=1 



= AA 3 (t a ,t 5 ) |l) 3 ^( 3 )(x 6 -x a -ptv(t)) 



• exp < z 



dt 



^v 2 (t) - V(x(t)) 



(4) 



Here the "measure" is given by 2? 3 -u = flj^ ^ 3y fc an d the normalization factor 



jV(t a ,t b ) := (^Jvvexp ij^dt^v 2 {t) 



(5) 



ensures that the Gaussian integral gives unity as is evident from the discrete form. Note that the 
functional integral over v does not require any boundary conditions which are all contained in the 
remaining 5 function. A more symmetrical form for the argument of the potential is obtained by 
writing 

x a + x 6 1 /■*(- 



x(t) = ^±^+ 1 -J t b dt' Sg n(t-t') V (t') 



(6) 



where sgn(x) = x/\x\ is the sign function. We have x(i) = v(t) and the boundary conditions for the 
paths are fulfilled due to the S function in Eq. (4). We now write Eq. (2) as 



Si^ f = lim e ^ Ei+E ^ T [ d 3 xd 3 ye- i ^^U(x,T;y,-T)e i ^ y 

T— »oo J 



(7) 



and insert the representation (4). Using the coordinates r = (x + y)/2 , s = x — y we then obtain 

r+T 



^ 2 (t)-K.v W 



Si^ f = Im^e^+ZfV J d 3 re-^ r M 3 (T,-T) J V 3 v exp | i j_ 

■exp | -i J ^ rftV(r + x„(t))| , 



(8) 



since the relative coordinate s is fixed by the 5 function in Eq. (4). Here we have defined the 
momentum transfer and the mean momentum by 



kf-ki, K = -(ki + kf) 



(9) 
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Furthermore, 



Mt) = ^J^dt'sgn(t-t')v(t') 



(10) 



where the subscript denotes the dependence on the variable over which one is integrating functionally 
[24]. The shift v(t) — ► v(t) + K/m eliminates the linear term in the exponent of the functional 
integral (8). Since 

J dt' sgn(i - t') = 2t for t £ [-T, +T] 

and 



Ei + Ej 



m 



Am 



we obtain 



. (ki-k/) 2 






4m 




Am 


J V 3 v exp 


'J 


f +T , m n, . 









(11) 

(12) 



x exp 



(13) 



Note that Eq. (12) is valid without energy conservation which has not been imposed (derived) yet. 
With no interaction we obtain 



:(0) 



i->f rp 

and therefore we will consider 



lim exp f i^-T] (2vr) 3 5 (3) (q) = (2tt) 3 <5 (3) (kj - k/) (14) 
»oo \ 4m y 



(5-1)^ = ^exp yi^T j y d 3 re- jq ' r AT 3 (T,-T) J V 3 v exp i J ^ dt^v 2 (t) 

i j + dtV (^r + ^t + x v (t) 



x < exp 



(15) 



in the following. Since the potential vanishes at infinity Eq. (15) is a well-defined integral. 



III. Asymptotic times: Elimination of dangerous phases 



The path integral representation (15) is exact but suffers in the present formulation from the explicit 
appearance of a "dangerous phase" q 2 T/(4m) proportional to T, in the first exponential of Eq. (15). 
It can be checked, of course, that this phase cancels in each order of perturbation theory so that the 
limit T — ► oo can indeed be performed but one would like to have a formulation where this phase does 
not appear at all. This can be achieved by recognizing that each power of q 2 arises from applying the 
the three-dimensional Laplacian —A to the factor exp(— iq • r) in the integral over r. An integration 
by parts then lets it act on the potential term [25]. In order to reduce it to a shift operator one may 
"undo the square" , for example by a 3-dimcnsional path integral 



exp 



4m 



-TA 



M* 3 (T,-T) [v 3 wexp -i [ +T dt — w\t)± [ +T dt-f(t)w(t)-V 
J J-T 2 J-T 2 



(16) 
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Here f(t) should fulfill 



J-T 



dtf(t) = 2T 



(17) 



and the normalization again ensures that the pure Gaussian integral gives unity. 

Note that the sign of the quadratic term in the exponent necessarily is reversed if one wants to 
have a real shift operator whereas the linear term can have any sign. Real arguments of the potential 
are mandatory if an analytic continuation of the potential into the complex plane is to be avoided. 
Such a procedure would depend on the specific analytic properties of the potential and would have 
to be considered on a case-by-case basis. We will call w(i) an "antivelocity" and choose the negative 
sign in the linear term for convenience. Also we will take 



/(*) = sgn(-t) 

so that the shift operator simply becomes 

^J^dtsgn(-t)w(t)-V 



exp 



= exp [-x U) (0) • V; 



(18) 



(19) 



and the antivelocity degrees of freedom are as close to the velocity ones as possible. Then we obtain 
the following path-integral representation for the S matrix 

(S - l) w = lim J d 3 r e- iq ' r \M(T, -T)| 6 J V 3 vV 3 w exp i J +T dt y (v 2 (t) - w 2 (t 



exp 



•; / dtV r + — t + x„(t) - x^fO) 
J-T V rn 



1 . (20) 



There is an interesting analogy with the Lee- Wick approach to Quantum Electrodynamics where also 
fields with a wrong-sign kinetic term are introduced [26] to remove all infinities. These "phantom" 
degrees of freedom are often described in an indefinite inner product space [27]. In our case, however, 
they are not conjectured but necessarily appear when eliminating the asymptotically diverging phase 
q 2 T/(4m) in the S matrix 

At first sight the present approach to remove the infinite phase looks as if the phase space path 
integral used in Ref. [18] has come back through the backdoor in disguise of a functional integration 
over velocities and antivelocities. However, since the argument of the potential in Eq. (20) only 
depends on the fixed quantity x lu (0) , the path integral over the antivelocity is not a full functional 
integral but could be replaced by an ordinary one. Such a representation corresponds to using 



exp 



T A 

4m 



J d 3 w exp (—imTw' 



d w exp 



-imTw 2 ±Tw V 



(21) 



i.e. a constant antivelocity w, instead of Eq. (16). This may offer definite advantages in all cases 
where an additional functional integration would be costly as in attempts to evaluate the real-time 
path integral numerically. However, compared to Eq. (16) it has the disadvantage that an explicit 
dependence on the time T formally remains and that v, w are treated differently. Therefore we will 
use the time-dependent antivelocity w(t) in the following applications. 
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IV. Faddeev-Popov methods for the T matrix 



How to extract the T matrix from the S matrix ? For weak interaction one can develop in powers of 
the potential and one finds that in each order an energy-conserving 5 function can be factored out. 

To achieve this without a perturbative expansion of the S matrix one can use the trick which 
Faddeev and Popov (FP) have introduced in field theory for the quantization of non-Abelian gauge 
theories as was first proposed in Ref. [18]: We note that in the limit T — > oo the action in the path 
integral (20) is invariant under the transformation 



t + T, 



K 



m 



v(t) = v(f) 



(22) 



since 



dtV (^r+^t + x v (t)-x w (ofj = T dtvft + ^+\ J d?v(F)sgn(t-F) -x„,(0) j . 

V -T-T ' 

(23) 

For finite r and T — > oo one may expect that the change in the integral limits is of no relevance 
and therefore that the action remains invariant under the above transformation. Actually the limit 
T — ► oo is nontrivial and needs a more rigorous investigation which is beyond the scope of the present 
investigation. Instead we will verify that our procedure is correct by checking that each term of the 
Born series emerges from our path integral representations. 

If the action is assumed to be invariant under the transformation (22) then it does not depend on 
the component of the vector r which is parallel to K, leading to a singularity when integrating over 
that component. This singularity is just the energy-conserving 5 function we are looking for. We can 
extract it by first fixing it and then integrating over all possible values: For example, we multiply the 
path integral (15) by the following factor 



i = M 


r+o° / 


K " 






/ dr 5 K ■ 


r H r 


m . 


' — oo \ 


m 





(24) 



where A is an arbitrary fixed ("gauge") parameter and K = K/|K| the unit vector in the K direction. 
We now perform the transformation (22) in the path integral and obtain 

IKI f +oc f ( K \ / * \ 

(5-l) w = ^ T lim dry^rexp^-iq.r + iq.-rj^K.r + A) 



\M(T,-T)f J V 3 vV 3 w exp | i ' dtj [v 2 (t) - w 2 (t) 



x < exp 



J + t dtV (v+ ^-t + x v (t) -^(Ojj 



1 



(25) 



To simplify the nomenclature the original variables are used again. The only dependence on r in the 
integrand now resides in the factor exp(— irq • K/m) and thus the integration over it produces the 
energy-conserving 5 function [28] 



2tt5 



q K 



m 



( k f k 2 s 
\ 2m 2m 



(26) 
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In addition, after the transformation the longitudinal component of r is set to the value —A. Noting 
that <7|| = we then obtain the following expression for the T matrix 



+00 

7^-3) = £ j d 2 be-^ h |AA| 6 J V 3 vV 3 w expj i J dt^ [v 2 (t) - w 2 (t)} 



—00 



o*XK(b,' 



(27) 

Here we have taken the limit T — > 00 and have written the corresponding Gaussian normalization 
factor as 

M := A/"(+oo,-oo) • (28) 

In Eq. (27) the phase xk is defined as 

/+00 / K - \ r+00 

d tv(h + -t + x„(t) - Xlu (0) - AK J =: — / <ft K(£ K (t)) (29) 

while b = rj_ denotes the transverse component of the vector r (the impact parameter). With 9 being 
the scattering angle, we have 



g=|q| = 2A;sin^-J , K = |K| = ^cos^-J . (30) 

Writing A = Kto/m we see that the "gauge parameter" can be traded for an arbitrary time to m t ne 
reference path b + K(t — to)/m . We expect that A = 0, i.e., to = is the most symmetric choice (see 
below) . 

As an exact path-integral representation of the T matrix Eq. (27) is one of the major results 
of this paper. The superscript "3-3" indicates that in addition to the three-dimensional velocity 
variable a three-dimensional antivelocity is used to cancel divergent phases in the limit of asymptotic 
times. Using Eq. (4) backwards it is also possible to write the result as an ordinary path integral 
over paths x(i),y(t) instead of velocities v(t), w(t). These paths have to fulfill boundary conditions 
x(±T) = ±xo/2,y(±T) = ±yo/2 and one has to integrate over xo,yo at the end. However, this 
brings neither simplifications nor new insights and so we will not pursue it further. Instead we will 
show in the next section that one can obtain the desired cancellation of divergent phases with an 
one-dimensional (longitudinal) antivelocity only. 

V. Ray representation 

The representation (27) can be simplified by a simultaneous shift of the impact parameter and the 
velocities 

v(t) = ^sgn(t)+v'(t) , w (t) = ^-sgn(t)+w'(t) (31) 
b = b'-x^ ± (0) + x^ ± (0). (32) 

This transformation is suggested by a stationary phase approximation to Eq. (27) 

5 r+ T 



5v(s) J-T 



r. 



dt 



fv 2 (t)-F(£ K (t)) 



, (33) 
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which gives for the stationary values of velocity and impact parameter 



mv stat (s) = J dtVV(£ K (t)) ^sgn(t-s) 
mw stat ( S ) = dtVV (H K (t)) ^sgn(-s) , 



<ftv 6 v(£ K (*)) 



(34) 



We thus find w^ at (s) = q/(2m)sgn(s) and for small scattering times [29] t or asymptotic external 
times s also 

vi tat ( S ) « ^sgn( S ) (35) 
which suggests the shift (31). However, doing so introduces additional terms in the exponent since 

= f dt v' 2 (t) - q • x„. ± (0) + £-T . (36) 



Similarly 



f /_ t ^w 2 (t) = |^ ^ w '2 (t) _ q . Xw , ±(0) + |Lr 



(37) 



so that 



+oo 
^ J — oo 



v 2 (t) -w 2 (t) 



+0O 



v /2 (t)-w /2 (t) +q-[x U) ,_ L (0)-x t ,_L(0)] (38) 



is independent of the time T used for regularization. But finite terms remain which are then canceled 
by the shift (32) of the impact parameter. Note that only the transverse component of x v > and x ro / 
can appear in Eq. (32) since the impact parameter necessarily is a two-dimensional vector. This 
asymmetry between perpendicular and parallel components can be traced back to the constraint (24) 
and will persist in the following formulae. Using the relation [30] 



f+T 



J dssgn(s-t)sgn(s-t') = 2 T-\t-t'\ , t,t' € [~T,+T] 



we find from Eqs. (10) and (31) that 



q 

2m 

q 



|*| -T +x„,(t), 



*v,(0) = ^(-T)+x^(0) 
Therefore the argument of the potential term also becomes (formally) T independent 

€k(*) - W*) = b' + ^-AK + x,(t)- Vl (0)-x,„,,|(0)K. 

Here 



m 



Pray (t) = K + |s g n(t) = k i e(-t) + k / e(t) 



(39) 

(40) 
(41) 

(42) 
(43) 



is the new momentum along which the particle mainly travels: for t < it is the initial momentum 
and for t > it is the final momentum. This is also what one expects intuitively at high energies 
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Figure 1: Scattering geometry for a potential of radius R, the impact parameter b, the ray made by 
the incoming and outgoing momenta kj f , and the mean momentum K = (kj + kj)/2. 



and is depicted in Fig. 1. Note that the magnitude of p ray (i) is k for all t and therefore the velocity 
of the high-energy particle along the "rays" remains the asymptotic k/m instead of the unnatural 
K/m = k cos(9/2)/m. 

After shifting of arguments the integrand does not depend on x w ±(0), i.e., wj_(t) anymore. There- 
fore the integration over the perpendicular components of w(t) can be performed trivially cancelling 
the corresponding Gaussian normalization constants. Choosing A = 0, omitting the prime for the 
shifted variables and writing w for wn the new path integral representation now reads 

T^/ ] =i^J d 2 be- l * h N z N* J V 3 vVwexp li J dtj [v 2 (t) - w 2 {t)] j j e W(b,v,«0_i j (44) 

with only one (longitudinal) antivelocity which is indicated by the superscript "3-1". The phase is 
given by 

+ 00 / v 

Xray(b,v,«;) = - J dt V I b + + x v (t) - X, ± (0) - x w (0) J . (45) 

Note that both path integral representations of the T matrix are not impact parameter representations 
in the strict sense since both the phases x an d the factor K = k cos (6/2) carry an angle dependence 
whereas in an exact impact parameter representation of the T matrix this dependence would only 
reside in the factor exp(— iq • b) [31]. As a consequence, unitarity of the S matrix, i.e., validity of the 
optical theorem is not immediately evident although these are exact path integral representations. 
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A. Microreversibility 

It is worthwhile to explore how microreversibility (time reversal) of the T matrix [32] is realized in 
the present path integral approach. This is the invariance under the exchange 

k . ^ _ k/ , k/ ^ _ ki ) ( 46 ) 

i. e., 

q^q , K^-K. (47) 

We first note that the gauge parameter A has to vanish since it multiplies the odd vector K in the 
argument £k(*) °f the phase 

A = . (48) 

This is also evident from the FP constraint (24), where the argument of the 5 function would have 
different parity upon time-reversal or simply by considering A as an arbitrary time scale to f° r the 
longitudinal motion which would destroy the time-symmetry between initial and final states. 

However, microreversibility does not constrain the dynamical variable v(i). Let us discuss that for 
the case of an one-dimensional antivelocity with the phase Xray(b, v, w) given in Eq. (45): the impact 
parameter b is unaffected, but the reference path obviously changes under the transformations (47): 

Pl^lt ^ --t+^-\t\. (49) 
m m 2m 

This can be compensated [33] by changing the integration variable t — > —t: 

Xray (b,v,u,) - - [ + °° dtV fe ray (-t)) (50) 

where 

£ ray (-t) = b + -t + ^-|t|+x„(-t)-x„ ± (0)+x„,(0)K (51) 

and 

1 r+oo 1 r+co 

M-t) = o dt ' ~ v (*') = o dt ' ^ - *') H v ("0 • ( 52 ) 

£ J -co £ J -co 

Decomposing the variable v(t) into even and odd components 

v(i) = v+(i)+v_(t) with v±(-i) = ±v±(i) (53) 
one sees that the kinetic term is quadratic in both components 

f+oo ffy. r+co 



y / °°dtv*{t) = ^ / °° dt U(t)+vl(t)} . (54) 

^ J —CO ^ J —CO 



This allows us to transform 

\ f+OO 



1 r+oo 

M-t) = ~ \ dt' sgn(i - f) (-) [v+(f)-v_(f)] • (55) 

& J — CO 

into x„(f) by a simple change of integration variables 

v+(t) - -v+(t) (56) 
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in the velocity path integral (leaving v_ unchanged) and demonstrates invariance of the phase Xray 
and of the whole T matrix. Of course, the subtraction terms 

»+oo ^ r+oo 

2 J 

depend only on the time-odd components. For the case of a three-dimensional antivelocity the argu- 
ments are even simpler but completely analogous. 



\ r+oo 

*v±(0) = / dt'sgn(-t')v ± (t') 

£ J — oo 



d£ sgn(-i') w(t') 



(57) 



B. Tests 

As a test for the correct treatment of the various limits and shifts which we have performed, the Born 
series should be obtained from the path integral representations T^ 3 " 3 ) and T^ 3-1 ) . Here we only 
consider the first Born approximation while terms of arbitrary order are evaluated in appendix A . 
The first-order T matrix is simply obtained by expanding the corresponding phase to linear order and 
Fourier transforming the potential 



(3-3) Born 



! ^ J d 2 be-^ b J ^3 V(p) \Nf j V 3 vV 3 w exp j J dt [v 2 (t) - w 2 (t) 



+oo , 

J ds exp I ip 



b + — s + x„(s) - x w (0) 
m 



The functional integrals here are simple Gaussian ones of the form 



G {d) := N d J V d v exp J i J 



dt 



T 



V(t) + g(t)-v(t) 



exp 



T 2m 



(58) 



(59) 



and we let the time T go to infinity only at the end of the calculation. From the relation (10) we read 
off g v (t) = psgn(s — t)/2 for the v integration and g w (t) = psgn(— t)/2 for the w integration. Thus 



+T 



sgn 2 (s -t)- sgn 2 (-t) 



exp -i|-(2T-2r) 



= 1 

(60) 



and 



T, 



(3-3) Born 



m J 

I 



d 2 be 



-iq-b 



f d 3 p ~ f+ T ( K 

/ -. — r-r V(p) lim / ds exp ip • b + ip • — s 

J (2vr) 3 t^oo L ^ V 



m 



d 3 P y(p)<5( 2 ) (p ± -q) -tff-j 



= V(q,0) = y(q) 



(61) 



as expected. 

Although the ray representation (44) was derived by a simple shift of integration variables from 
Eq. (27) and therefore did not involve any additional large-T limits it is instructive to derive the first 
Born approximation explicitly in this case too. We have 



r (3-l)Born ; 



K 



m 



+0O 



Pray(^) 



+oo , 

x J ds exp< ip • 



b + 



m 



dt v 2 (t) -w 2 (t) 
s + x 1) (s) - x t) _ L (0) - x w (0) K 



(62) 



-oo 
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From the master path integral (59) we obtain 

G^G^* = expj-^- [pi( S gn( S -t)- S gn(-t)) 2 +pj( S gn 2 ( S -t)- S gn 2 (-t) 

(63) 

corresponding to perpendicular and parallel path integration over the velocities v, w. The longitudinal 
component of the momentum p is completely cancelled by the contribution from the antivelocity w 
but now a term remains in the exponent which is proportional to p^. Performing the t integration 
by means of Eq. (39) all T dependence cancels and we obtain 

G (<*=3). G (d=D* = exp{-^-pi|,|} . (64) 
Using the explicit form (43) of the momentum p ray (s) it follows that 

-7ft 1)B ° rn = - U%e^ [-^V(P) 
t_>/ m J J (2tt) 3 

/+oo / K q i \ 
ds exp ( ip ■ b + ip ■ — s + ip ■ \s\ p 2 [ \s\ j . (65) 
-oo V m 2m 2m ) 

The b integration leads to pj_ = q and therefore the leftover term from the vj_ integration is taken 
away by the contribution from the modified reference path. Thus we obtain again the correct first- 
order result (61). 

In Appendix A we show how to obtain the complete Born series from these two path integral repre- 
sentations. This demonstrates that they are completely equivalent to the standard (time-independent) 
scattering theory and can be utilized without doubt. 



VI. High-energy expansions 

The path integral representations (27) and (44) for the T matrix are the natural starting points 
for high-energy approximations. Under these kinematical conditions one expects that the particle 
essentially moves along straight lines with a constant velocity and that the functional integral over 
velocity and anti-velocity only describes the fluctuations around this trajectory. 



A. Eikonal expansion 

Taking Eq. (27) (where the particle travels along the mean momentum K) as reference one indeed 
finds that this is the case: By setting 



m 



t = - z , v(t) 



m 



v(z), w(t) 



m 



w(z) 



it is seen that the path integral (27) takes the form 

.K 



{+oo 
% - J dz [v 2 (z) - w 2 (z) 



exp 



+oo 

m f 

~ 1 kS c 



dzVl b + Kz + 



xs(z) - Xtf(O) 



1 . 



(66) 



(67) 
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In many applications (e. g., in atomic physics) the energy of the incoming particle is not large 
compared to its rest mass. Therefore we consider m/K not as small but as fixed in the following. 
Equation (67) shows that this factor just multiplies the potential but - irrespective of its magnitude - 
a systematic expansion in inverse powers of K of the T matrix is possible for fixed momentum 
transfer (or scattering angle). This is achieved just by expanding the phase simultaneously in powers 
of v(t),w(£) and performing the functional integral term by term: At high energy the fluctuations 
around the straight-line trajectory are indeed small. Of course, the convergence will depend on size 
and smoothness of the potential as higher and higher derivatives of it will appear in the expansion. In 
addition, since K = k cos(#/2) becomes smaller in backward direction the convergence of the expansion 
will deteriorate for larger scattering angles. A rough estimate of the validity of the expansion may be 
given by the requirement that the next order term of the Taylor expansion be small compared to the 
leading term 

<C V . (68) 



! W • x c 



'K 

Assuming that the velocity fluctuations are only relevant within the range R of the potential one finds 
v = 0(1/ y/R) and x, d = 0(y/R) and thus 

where a is the scale over which the potential changes appreciably. 

Let us start with the lowest order term. Setting v = w = in the argument £k(*) °f the potential 
immediately gives 

i+oo 
' X%0) = -%JdzV(b + Kz), (70) 
— oo 

because the functional integrals are trivially one by normalization [34] . This is a variant of the eikonal 
approximation due to Abarbanel and Itzykson [35] where K = fccos(#/2) appears everywhere instead 
of the asymptotic momentum k. For a spherically symmetric potential V{r) we have the standard 
result 

X ( >) = ~ j™dzV(VvT*) • (71) 

It is easy to calculate the next-to-leading order correction by expanding the phase up to linear 
order in v(i) and w(i), and performing the shifted Gaussian integral by means of Eq. (59). The result 
is 

T^f = T« = <| J d% e~^ h | exp [i x <g + i x <»] - 1 j (72) 

with an additional phase function 

1 f 

xS(b) = -— lim / ds / dhdt 2 Wi • VV 2 sgn(h - s)sgn(t 2 - a) - sgn 2 (-s) 

1 r+oo 

= / dhdt 2 VFi • W 2 \h-t 2 \ (73) 

4?7i J — oo 
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where Wi is an abbreviation for W(h + Kzj). Again the contribution from the antivelocity w 
naturally cancels explicit T terms when the integration over s is performed with the help of Eq. (39). 
In appendix B it is shown that for a spherically symmetric potential the expression simplifies to 



x>) - ~ If 



1 



m 



2 r 



1 



b db 



f 
Jo 



dz V 2 ( 



(74) 



This is identical with the phase T\{b) in the systematic eikonal expansion of Wallace [14] apart from 
the appearance of K = kcos(6/2) instead of k which is unimportant in this order and for forward 
direction. Note that this additional phase already appears in exponentiated form as conjectured by 
Wallace. 

One may wonder whether Eq. (74) is the correct result up to order K~ l because the next order 
term is also of that order. For insight into this question it is instructive to consider the example of 
an one-dimensional (ordinary) integral 

/ dx exp(ix 2 ) exp [— iV(a + \fex)\ 



T(a,e) 



= / e iV(a+^x) 



(75) 



/ dx exp(ix 2 ) 

After expanding the function V in the exponent for small e as V(a + yfex) = V(a) + yfeV'{a)x + . . . , 
keeping terms up to order x 2 in the exponent, expanding higher-order terms, integrating term by term 
and re-exponentiating one obtains 

T(a,e) = M J dxexp\ i -iV(a) + ix 2 (l-eV"(a)/2)-ix^V'(a)\ • [l + 0(e 3/ V)] 



exp 
exp 



-iV(a) - 
-iV(a) 



IE 



A{l-eV"{a)/2) 



V' 2 {a)--Hl-eV"{a)/2) 



■{l + 0(e 2 )} 



l lv' 2 {a) + - A V"(a) + 0{e 2 ) 



(76) 



Since the correction phase Xai was obtained by truncating the Taylor expansion of the potential at first 
order it just corresponds to the second term and we seem to have missed another, purely imaginary 
phase linear in the potential which is also first order in e or 1/K . However, closer examination shows 
that this is not the case. It is, of course, possible to prove that assertion directly by evaluating the 
required functional integrals. These are more general Gaussian integrals of the type 



M J Vv exp I i J dt 



\t) +y/lg{t)-v{t) 



+ e 



j dtdt' v(t)h(t,t')v(t') 



exp 



(77) 



and multidimensional extensions thereof. However, there is an easier approach using the cumulant 
expansion (see, for example, Ref. [36]) which in the one-dimensional example of Eq. (75) reads 

„-2 



T(a,e) 



exp 



ici(a,e) + c 2 (a,e) + . . . 



ci(a, e) = e ( V(a + yfex) ) 

c 2 (a,e) = e 2 ( (V(a + Vex) - (V(a + yfex))' 



(78) 

(79) 
(80) 
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Of course, by expanding the cumulants in powers of e one obtains the same result (76) as before. 
Application to the eikonal expansion is straightforward: it is easy to calculate the cumulants in closed 
form and since Xaj 1S quadratic in the potential we only have to expand the first cumulant in inverse 
powers of K in order to obtain all terms which are linear in the potential. This is very similar to 
working out the first Born approximation and we obtain 



d 3 p 
(2^)3 

d 3 p 



+00 



V(P) 



ds exp 



. K 

ip • I b H s 



+00 



ds exp 



ip 



(2k) 

%jy°dzV(b + Kz) EE xssw . 



m 

. K 

b H s 

m 



exp [ip • (x^(s) - -x w (0))] 



(81) 



where the (functional) average over v,w with the weight exp[im / dt (v 2 — w 2 )/2] gives one according 
to Eq. (60). Hence there are no higher-order terms linear in the potential beyond the leading eikonal 
phase and Eqs. (72) and (73) are correct up to and including order 1/K. 



B. Ray expansion 

The path integral representation (44) gives rise to a different high-energy expansion because we expand 
around the momentum p ra y(^) which takes into account the different asymptotic directions before 
and after the scattering. While this complicates the analysis and leads to an additional momentum 
transfer dependence some advantages at larger scattering angles may be expected: Applying a similar 
scaling argument as in Eqs. (66, 67) one sees that now a systematic expansion in inverse powers 
of k is obtained which we will call the "ray" expansion. A disadvantage is the expansion around 
a discontinous reference path which abruptly changes direction at t = 0. This may deteriorate the 
convergence properties of the expansion but may be remedied by another choice of the function f(t) 
in Eq. (18) subject to the constraint (17). 

The lowest order term is obtained by setting v = in the argument of V and immediately gives 
a new high-energy approximation 



iqb 



• (0) 
g*Xray 



(82) 



with a phase [37] 



X^(b,q) 



m 
T 



+00 

J dzV(p(z)), p{z) = b + 



p(t = mz/k) 



k 



v. q , K 



-00 

roo 



ee -jj^dz [V (b-M) +v(b + i f z) 



(83) 



This has some similarity with the eikonal phase derived by Levy and Sucher [17] although these 
"symmetric" eikonal expansions [38] are quite different from our approach. 

For a spherically symmetric potential we have for the leading order ray phase function 

X^ y (bJ) = 



dz V(p(z)) 



(84) 
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where 

p(z) = ^/b 2 + z 2 + h-q_z/k = ^b 2 + z 2 + 2bz(3 . (85) 

Here we have defined 

P = ^ = srnf 6 -) costp, |/3|<1, (86) 



2k \2, 

where ip is the angle between the impact parameter and the momentum transfer. In forward direction 
(where all different eikonal approximations should be equivalent) this is seen to reduce to the usual 
eikonal phase plus a correction: 

*ffiM>^-£f*V W -£^f *>) + .... (ST, 

With dV(r)/dz = zV'(r)/r the correction term can be easily integrated and gives 

Xffi(b,/J) ^ - 2 -^^ dzV (r) + h.^V(b) + ... . 

Combining the result with the exp(— iq- b)-factor in the impact parameter integral it is thus seen that 
the main effect of the ray approximation is the replacement of the momentum transfer by an effective 
momentum transfer 

q e s(b) = q(l-^V(b) 

which takes into account the energy gained (or lost) by moving in the attractive (or repulsive) potential 
at closest approach: 

k 2 ^ k 2 s (b) 



2m 2m 



+ V(b, z = 0), q cS = 2k cS sin (0 . (90) 



This approximation (with an average, constant value of the potential) is standard practice in electron 
scattering from nuclei where higher order effects are roughly included by evaluating the Born approx- 
imation form factor as function of an effective momentum transfer [39]. However, when doing that it 
is also well known [40] that a flux factor (k e g/k) 2 is needed for the scattering amplitude. 

This flux factor is provided by a purely imaginary phase Wriy which appears in next-to-leading 
order and corresponds to the second term in the example (76). In contrast to the eikonal expansion 
in the previous subsection this correction does not vanish anymore. Let us evaluate it by calculating 
the first cumulant: 



(2vrj3 V ^ J-oo dS 6XP ^~ iP ' p(yS ^ 



(2vr) E 

x ^exp|-ip- [x v (s) -x„_l(0)] +ip-Kx w (0)| ^ . (91) 
The average in the last line has already been evaluated in Eq. (64) so that 

[ p 2 

— \f(n\ pyti f — in • nf?.}] pyti — i- 

(27T) 

To order k~ l there is now a purely imaginary phase with magnitude 

^r ( S(b,q) = / J^^^pl ex P[-^P-pW] = -^kJ A <> dz\z\V(p(z)) , 

(93) 



(x ray ) = -^/ + J^|^3y(p)exp[-,p.p(z)]exp -i^L\z\ =: *g +0 (AT 2 ) . (92) 
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where 

b ~ db 2 + bdb + b 2 dip 2 bdb db + b 2 dip 2 ' ( ' 

For a spherically symmetric potential this simplifies to 

1 Til Z* 00 

uW(b,0,0) = -- — A b J Q dzzV(p(z)), (95) 

Note that Wray now depends on three variables (apart from the overall powers of 1/k): b, b • q, q or 
b, (3,6 . This is because the Laplacian A& contains explicit derivatives with respect to ip . 

What is the effect of the real factor 

e-"™y ~ 1 - (96) 

on the scattering amplitude? If the improvement from the leading ray approximation is incorporated 
into an effective momentum transfer (as discussed above) we may consider exp(— 

J&) simply as an 

amplitude correction as discussed above. However, an alternative interpretation arises if the correction 
(88) is included in a scaled impact parameter 

* = »(l-£v W ) => » = y(l + ^V(«0) + o(^;). P7) 
This implies the following change in the integration measure 

bdb = b'db' 1 + p ( F(ft') + A (6V(6') )) • (98) 

However, in forward direction p 2 (z) r 2 = b 2 + z 2 and therefore 

m m 1 d , f°° , dV m 1 d , f°° , , dV m r T ,, IN , T , //Ixn , N 



w ray 



so that exp(— w^y) exactly cancels (at least in the forward direction) the Jacobian arising from the 
scaling transformation. The leading order ray phase with the scaled impact parameter as argument 
is 

XS(O-O) * - 2 ^fdzV(r')-^V(b')fdz^l 

/-oo 2m 2 (9 Z" 00 , 

= / dzV(r')-^V(b')b'^- dzV(r'), r' = VV 2 T^ 2 (100) 

which has a correction term similar to Eq. (74) in the eikonal expansion. Thus the leading order ray 
expansion already contains approximately higher order eikonal terms. 

Of course, there is also a real first-order phase Xray which is obtained by expanding Xray up to 
first order in v, w : 

Xr, y (h, V ,w) = X f2?+ j^dt [g v (t)-v(t)-g w (t)w(t)} + ... (101) 
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where now 



[&,(*)]* = "J /_ T WPi) [ s S n (*i "*)-(!- M sgn(-t) ] 
1 /"+ T 

9w{t) = -J dti a fc y(p 1 )sgn(-t) . 



(102) 
(103) 



Here fc = 1,2,3 are the cartesian coordinates of the vector and the argument of the potential is 
always p 1 = b + x re f(ii) = b + Kti/m + q|ti|/(2m) . Applying the Gaussian integration formula (59) 
we obtain the real correction phase of order one 



1 

2m 



+T 



dt 



ll{t)-g 2 w {t) 



(104) 



and performing the t integration with the help of Eq. (39) we find - as expected - that all T 
dependence cancels. Thus 

xS(b,q) = -^J^dhdt? {W(*i)-W(t 2 ) \ti-t 2 \-V h V(h)-V h V(t 2 ) (\h\ + \t 2 \)} . (105) 

For a spherically symmetric potential some algebra which is outlined in appendix B leads to 



1 m 2 

kip 



/ f) roc 



(106) 



where p and (3 are defined in Eq. (85) and (86), respectively. The first term in Eq. (106) is identical 
with Wallace's eikonal phase t\ for forward scattering when p — > r. It may be surprising that Eq. 
(106) contains an additional term which does not vanish in forward direction, i.e., for (3 = 0. But 
this is just the term which exactly cancels the last term in the approximation of Eq. (100) so that 
the correct first-order eikonal expression for forward scattering is obtained. Note that there is no 
singularity in Eq. (106) at (3 = ±1 as can be also seen in appendix C. We therefore have in first-order 
ray expansion 



-n-*/ — -'ray 



- / d 2 be~^ h f 
m J L 



exp 



jy(0) +iy W _ w (l) 

"Aray 1 "Aray ray 



1 



(107) 



VII. Numerical results 

Let us test the high-energy expansions for the case of scattering from a Gaussian Potential 

V(r) = V e- r2 / R2 



(108) 



with the parameter values 2mVoR 2 = —4, kR = 4 , i.e., E = — 4Vq , corresponding to the case 
where a particles scatter elastically from a particles at 166 MeV center-of-mass energy (R = 1 fm). 
The parameters are precisely those where convergence of the standard eikonal expansion was found 
to be unsatisfactory [14]. For completeness the analytical expressions for the various phases of a 
Gaussian potential are listed in appendix C. We have evaluated Tai and T ray by numerical integration 
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20 40 60 80 100 120 140 160 
6 (degrees) 

Figure 2: Differential cross section from a Gaussian potential with strength 2mVoR 2 = —4 at kR = 4 
as function of the scattering angle 9. Shown are the exact result from a partial wave calculation and 
from the zeroth and first order of the high-energy expansions derived in Sec. VI. In these the particle 
travels along the mean momentum [an eikonal approximation due to Arbabanel and Itzykson (AI)] 
or along a ray made up of the initial and the final momentum. 

using Gauss-Legendre rules with 72 points and a sufficient number of subdivisions of the integration 
interval which was mapped to a finite range by yi = R tan ipi where y = b,z for the AI expansion and 
y = b x , by, z for the ray expansion. 

Figure 2 shows the differential cross section obtained from these high-energy approximations compared 
to an exact partial wave calculation. The (AI) eikonal expansion shows the well-known failure at larger 
scattering angles and the corrections only slightly increase the point of deviation [41] . Since the cross 
section is sharply peaked in forward direction the total cross section is always well reproduced despite 
the deviations at higher scattering angles and not suited as a measure of (dis)agreement. 

The ray expansion does better at higher scattering angles at the price of being more complicated 
and less precise at small scattering angles. As scattering from a Gaussian potential at larger scattering 
angles is known to be dominated by many small scatterings these deficiencies may be attributed to 
the sudden change at t = which imparts a large momentum transfer to the scattered particle. In 
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addition, as was mentioned before a derivative expansion about this discontinous path will probably 
run into problems. It may be expected that a description based on a smooth path will do better but 
this will not be pursued in the present work. 

Figure 3 shows how the different expansions describe the cross section at higher energy (kR = 6). 
Again the ray expansion is closer to the exact result at higher scattering angles. 




20 40 60 80 100 120 
6 (degrees) 

Figure 3: Same as in Fig. 2 but for kR = 6. 



VIII. Summary and outlook 

Time-dependent methods for scattering have been investigated by several authors [42] in ordinary 
quantum mechanics. Using path integrals over velocities I have derived two new representations for the 
nonrelativistic T matrix which in a very natural way describe the propagation of high-energy particles 
in a local potential. Although the time evolution of the scattering process is also central in the present 
approach it leads to formulations which are quite different from the previous ones. This is because two 
important requirements must be fulfilled for obtaining a path integral formulation of the T matrix 
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from the matrix elements of the time-evolution operator U (T, — T) for infinite scattering times T. First 
one has to make sure that phases are eliminated which diverge for T — ► oo and, second, a suitable 
constraint has to be found which leads to energy conservation in the S matrix. In the present paper 
the first requirement is met by introducing phantom degrees of freedom ( "antivelocity" ) which cancel 
these divergences in a way reminiscent of the Lee- Wick proposal for Quantum Field Theory. Energy 
conservation is achieved by using the classic Faddeev-Popov procedure such that the component of the 
position vector parallel to the mean momentum K = (kj + k/)/2 is fixed. This involved some delicate 
(and at present not very well-controlled) limit procedures but we have checked that the resulting 
path-integral formulations of the T matrix produce the correct Born series in all orders. 

One of the advantages of these new path-integral formulations is that they can give rise to new 
approximation schemes or expansions. As they are close to a geometrical picture of scattering where 
the path integral describes the quantum fluctuations around some reference path it is not surprising 
that the eikonal approximation (in the variant of Abarbanel and Itzykson [35] where the particle 
travels along a straight-line path with velocity K/m) immediately follows and that corrections to it 
can be calculated systematically. A suitable scaling of variables in the path integral shows that these 
corrections involve inverse powers of K = kcos(6/2) and therefore inevitably grow at larger scattering 
angles 0. However, one is also naturally led to a new variant (ray approximation) which displays the 
different asymptotic directions along which the particle propagates at high energy and which should 
work at high energy irrespective of the scattering angle. Indeed, for high-energy scattering from a 
Gaussian potential some improvement over the Abarbanel-Itzykson eikonal expansion was achieved. 

There seems to be considerable room (and need) for improvement: a better Faddeev-Popov con- 
straint should eliminate the rather asymmetric treatment of longitudinal and perpendicular variables. 
It is unclear how Wallace's eikonal expansion [14] (where the particle is travelling along the mean 
momentum but with velocity k/m ) could emerge naturally from a path- integral representation. 
This formulation gives an impact-parameter representation of the T matrix and produces the exact 
Coulomb amplitude in lowest order [43] - which is not the case for the present formulation. A better 
control of the delicate limit T — > oo needed for obtaining the S matrix is certainly desired and finally 
one may ask whether a formulation without phantom degrees of freedom is possible. 

However, despite these shortcomings and the long list of desiderata our formulation seems to 
have some merits: at least it enlarges the "tool-box" of scattering theory and offers new possibilities. 
Among these one may expect new approximation schemes and, hopefully the prospect of evaluating the 
real-time path integral numerically, i.e., achieving a stochastic evaluation of the scattering process. 
Obviously this would be of great importance in the many-body case where one may assume the 
interaction potential as 

N 

V(r) = £V(r-r fe ) , (109) 

k=l 

with rk denoting the position of the feth scatterer. It is amusing that the path integral representations 
discussed in this paper lead to a multiple scattering expansion with exactly N terms when 

N N j 




n = i+E e n[ ex p(^)-i] (no) 



k=l j=l k 1 <k 2 <...kj 1=1 



is used. This is in contrast to Watson's multiple scattering expansion (see, e.g. ref. [44]) which 
contains infinite many terms and much closer to Glauber's theory where the incident particle cannot 
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scatter back due to its straight-line propagation. However, because of the subsequent path integration 
over all velocities the present expansion (if taken to full order) is not an approximation but allows for 
repeated scattering from the same scattering center. 

Extensions to relativistic scattering [45] also seem possible. Further investigations of this formu- 
lation as well as numerical studies of the real-time path integral will be reported elsewhere [46]. 
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Appendix A: Complete Born series from the path-integral repre- 
sentations 

Here we show that the various path-integral representations for the T matrix exactly reproduce the 
conventional Born series to all orders if the exponent is expanded in powers of the potential 



cxp 



+T 



dtvm)) 



n=l 



-iy 



YT 



dt 1 ...dt n V(Z(t 1 ))...V(Z(t n )) (Al) 



and the functional integrations are done term by term. This can be done by Fourier transforming the 
potential 

Vm)) = I 0f V ^ eiPvi{U) ( A2 ) 
We do that first for the version with a three-dimensional antivelocity as given in Eqs. (27), (29) where 
the reference path is along the average momentum 



K 



x^ffo) = —U and Z(U) = Zk(U) = b + x^)+x^)-x ro (0). 



We then obtain 



with 



T, 



■ £-4 



(A3) 
(A4) 



n=l 



K (-i) n 



n (/_>/ <gp(P. 



m n 



exp *X>-[b + x ref &)] G* 3_d) . 



(A5) 



i=i 
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For d = 3 we have to evaluate 



Gt 3) = \MfJ V 3 vV 3 w exp i J^dtj (v 2 - w 2 ) exp jzf> • [x^) - x. 



(0)] 



Since 



x^ti) - x w (0) 



1 f+ T 



dt 



sgn(ti - t)v(i) - sgn(-i)w(i) 



(A6) 



(A7) 



the path integrals to be evaluated are just Gaussian integrals of the same form as in Eq. (59) giving 
the result 



+T 



dt 



sgn(ti - t) sgn(tj - t) - sgn 2 (-t) 



(A8) 



Note that the first term in the square bracket comes from the functional integration over v and the 
second one from the functional integration over the antivelocity. As usual we have regulated the time 
integration by a finite time T which we finally will send to infinity. Using Eq. (39) any divergence in 
this limit is cancelled by the contribution from the antivelocity 



Gl 3 3 ) = exp { ~ £ ^ ■ Pj (T - \U - t 3l 



(A9) 



as was expected. 



Next we consider the ray representation (44) with an one-dimensional antivelocity which is a little 
bit more involved: first, the reference path is 



x ray (i) = Pray ^ = K f + _q_ 

m m 2m 



Second, the path integrals to be performed are again of the form (59) but with 

1 n 

gv±(t) = -^2p±i[sgn(ti-t)-sgn{-t)} 
z i=i 

■y n 1 Tt 

fl'w || (*) = T, 5ZP||*sgn(ti -t) , g w (t) = - ^p||jSgn(-t) . 



i=l 



i=l 



Therefore the Gaussian integration gives 
G { t 1] = exp 



r+T 



oW dt 

2m J-t 



S 2 v ±(t)+g 2 vll (t)-g 2 w (t) 



and after performing the t integral by means of Eq. (39) and some algebra one obtains 

= exp 



4^ [ p * ' p i I** ~~ *•?'! ~ P- 1 * ' P-Lj ( 1**1 + l*il 



(A10) 

(All) 
(A12) 

(A13) 



(A14) 



Compared to Eq. (A9) there is an additional term which, however, is exactly cancelled by the 
additional term from 



n /, \ n tv- n 

Pi 1 U = > P i ■ — U+ y,Pi ■ 7T~ ** 

m m r— f 2m 

i=i i=i i=i 



(A15) 
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if one takes into account that the b integration enforces 

n 



(A16) 



Thus in both cases the rath order term in the Born series reads 



x exp 



£Pi • ( b + -U) + ^ ^ ■ Pi l*i " *i 



4m .4 , 

i=i i,j=i 



(A17) 



where now the limit T — > oo has been taken. For further progress it is essential to recognize that 
the integrand is fully symmetric under exchange of t{ <-> ij as can be verified by the corresponding 
exchange p^ <-> Pj. If this is the case then 

n / r+oo \ r+oo /■<„ /-*2 

[J ( / d *» ) Fs y mm (*!■•■ = n! / rf i„ / dt n -i ... j dti F symm (h . . . t n ) . (A18) 

\J—oo / J —oo J — oo J —oo 

The factor in front cancels the factorial in the denominator of Eq. (A17). Furthermore, since the 
integration times are now ordered the last term in the exponential factor of Eq. (A17) becomes 



n n n n 

— Y, P* ■ Pi l*i - *il = ^ ^ Pi ' Pi (tj ~ ^ = 2^ S *J p J " 12 sgn ^' " fc ) Pfe ' ( A19 ) 

j,j=l i<j j=l k=l 



where sgn(O) = by convention. With the abbreviation 



1 



K 



:= ^ Pi " J2 s Z n ti ~ k )Pk + Pj ■ — 



k=i 



rn 



(A20) 



the time integrations can now be performed successively: 



f 2 dh e itl(ui - i0 ^ = 

J —oo 

f A dt 2 e it2{ui+U2 ~ i0) = 

J —oo 



'2 gi*2(ui-i0) 



u\ — iO 



2. git3(tti+«2— iO) 



u\ + U2 — iO 



/tn 
^ e it n - 1 (u 1 +U2+—+u n -i-iO) 



/+OO v-^n / 

Vi=i 



ui + u 2 + . . . + u n -i - iO 

n 



a it n (ui+U2 + ...+U n -l-iO) 



(A21) 



Note that the prescription how to handle the singularities arises from the requirement that the time 
integrations should converge at the lower limit. Alternatively one could give the particle mass an 
infinitesimal imaginary part m — ► m + iO already in the path integral so that exp (i J^t dtmv 2 /2) is 
damped. 
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Performing the b integration in Eq. (A17) gives another 5 function so that 



.K 



k=i 



d 3 p, 



n-1 



Uj 



(A22) 



u'=i 



Recalling the definition of Uj in Eq. (A20) we see that 



n n j£ n K 

2 n i = E Pi • — + ^ E s S n (i " fc )Pi • Pfc = E Pi • 



. m 2m . 

j=l 3=1 3,k=l 

because the last term changes sign under the exchange i <-> j. Thus 



m 



(A23) 



=1 V-E*=i«i + *o, 



P.7 



(A24) 



For n = 1 the last product is empty and the standard first-order Born approximation is obtained as 
already discussed in Sec. B.. 

For n > 1 the denominators in Eq. (A24) can be rewritten as 



k 1 ^ In 

Y. u 3 = ^z: Epj • E p*" 2K 

3=1 Zm 3=1 \i=k+l 



(A25) 



and the 5 functions allow us to replace 



E P«i = <1-Ep«-L ' E Pi 

i=k+l i=l i=k+l 



k 

E^ i 

i=i 



(A26) 



After some algebra one then obtains 

k 



E^ = ^ 



with E 1 = k 2 /(2m) = k?/(2m) = k^/(2m). This suggests the transformation of integration variables 
to 

Ik k \ 

4 := 5> l± -^,5>||+in , (A28) 



2m 



(x> ll+ *) 2 +(x>.,-fy 



(A27) 



i=l 



so that 



Pfe = 4 -4-i , k = 2...n. 
This also holds for k = 1 if we define 



2' 



(A29) 
(A30) 
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i.e., set empty sums to zero in the definition (A28). It is easily seen that the Jacobi determinant of 
this transformation is one. Furthermore the two 5 functions fix 



f 



K + - = k f 

2 / 



(A31) 



We then obtain the final result 
% = 



d 3 £ n -i ~ n ~ 2 

V {£-n f-n—lj 



(2vr) 3 

d 3 e n -i 

(2vr) 3 



n 

k=i 



d 3 h 



V(kf-£n-l 



(2vr) 3 E - Zl/ {2m) + iO 
1 



V(ti - to) 



(2tt) 3 
...■V(£ 2 -£i) 



E 



f 2 

2m ^ LU 



V (€ n _l - l n -2) 



ex 

2m 



V-(€i-ki) • 



(A32) 



This is identical with the usual quantum-mechanical expression obtained in time-independent scat- 
tering theory (here operators are denoted by a "hat" ) 



4> f v{x) 



V(x) . . . V(x) 



V(x) 0, 



(A33) 



E -p 2 /(2m) +i0 w "' y J E-p 2 /(2m)+i0 

when evaluated in momentum space [note the convention (A2) and the normalization of the scattering 
states 4>}. 



Appendix B: Evaluation of first-order corrections for a spherically 
symmetric potential 



For the eikonal expansion we start from Eq. (73), change to z = Kt/m and use VV(r) = rV'(r)/r 
where the prime indicates differentiation with respect to the argument. This gives 

4K3 L dz ^ 



n 



b + Z\Zi 



\z\ - z 2 \ ■ 



(Bl) 



Since r = \/b 2 + z 2 is invariant under z — > — z and the integrand is symmetric with respect to zi <-> Z2 
one obtains 



2m z 



CiZl 2i / CfZ2 — - — - 

r\ Jo V2 



The simple relations 



V'{r) _ dV(r) b V'(r) 



dV{r) 



(B2) 



(B3) 



r dz r db 

can be used to reduce the first-order eikonal phase to one-dimensional quadratures. The first one 
together with an integration by parts leads to 

,2 



2m z 
-£3 







dzV{r 



V'(r) 



b 2 -z< 



Then the second relation in Eq. (B3) may be employed to yield 



dz 



V(r)^ - zV(r) aV(r) 



dz 



(B4) 



(B5) 
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Finally another integration by parts in the last term gives Eq. (74). 



For the first-order ray correction the algebra is a little bit more involved. We start from Eq. (105) 
and use 

V'(p) 

u -r - — h -r 

m 



VV(p) = 



P 



(B6) 



Restricting the integration region to positive values one obtains 
1 f 



Pi 



P2 



,9 b • q . , k 2 

b 2 + -r- 1 (h + 1 2 ) + —2*1*2 



2m 



h + t 2 - \h-t 2 \ 



2m 2 



/ dZ\ / dZ2 Z2 + pb {Z\ + Z2) + Z\Z 2 

JO Pi JO P2 L 



V'{P2) 



fc d JO Pi JO P2 L J 

where in the last line the transformation z = kt/m and the symmetry of the integrand have been 
used. f3 is defined in Eq. (86). As in the eikonal case we would like to reduce this expression 
to one-dimensional integrals but the derivatives of the potential with respect to z, b are now more 
complicated: 

dv( P ) _ _ v'(p) , aui dv(p) _ _ v'(p) 



(Pb + z) , 



dz p db p 

Solving for zV'/p and bV'/p one now has instead of the relations (B3) 



(b + pz) 



V'(P) 
P 

V'(p) 



1 



As a final combination one needs 



1-P 2 
1 

T^p 2 

dV (5b 



dV 
dz 
dV 



P 



dV 
~db 
dV 



db 13 dz 



dV dV 
!h~ f ~db 



(B8) 

(B9) 
(BIO) 

(Bll) 



p dz 1 - p 2 

Equation (B7) multiplied by —k 3 /(2m 2 ) consists of four terms which can be simplified with the 
help of the above relations and appropriate integrations by part. The first term is 



a I J V'(pi) f, 

dz 2 Z2 



7 





dz\ 



V'{P2) 



Pi 



P2 



2\2 



oo 

j dZl 




P 



dV 

dz\ 

dV x 



dzi 



f dV 
\dz~2 



P 



dV 2 
db 



,8V 



dz\ 



dVo 



+P Z T— / dZ2 



db 



2{l-p 2 ) 2 
The second one reads 



2^ 



(B12) 



Jo 



pi 



I dZ2 Z2 

Jo 



P2 



Jo 2dz\ \Jo P2 

- V ^- l3 8il dzV ) 



1 pb_ 

2(1^ 



(B13) 
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Finally the third and fourth terms combined give 

V'(pi 



fib / dz x 

Jo pi Jo 

00 dV f Zl 



Z1 , 2 V (P2) , /X 
d,Z2 Z2 h 



P2 



dz\ z\ 



V'(p 







>V'( P 2) 



dz x ^— I dz 2 zi- 

OZi Jo P 2 



OC 



Pi 

V'(p) 



1) p 

Jo 



dz2 z\ 



V'{P2) 



pi 



dz Vz 2 

P 



dzV 2 . 



2 1 - (3 2 db jo 

Summing up all contributions we obtain the result given in Eq. (106). 



(B14) 



Appendix C: Eikonal and ray phases for a Gaussian potential 

Here we list the analytical expressions for the phases of a Gaussian potential 

1 



V[r) 



a = 



R 



For the AI eikonal expansion we have the well-known results [14] 
XAl[b) ~ K a 



For the ray expansion we get in zeroth order 



k a 



where 



F(B = abf3) := e B \iic{B) , F(0) = 1 



and erfc(x) = 1 — erf(x) is the complimentary error function [47]. The correction terms are: 



1 fmVo^ 2 



2k V k 



^1 e -2a 2 b 2 

2 a 



1 - Aa 2 b 2 (l - 



F 



+2v / 2a 2 6 2 (l ~/3 2 ) F(B) [2 - y/n BF(B)] j 



2 mVo b 2 
k k 



cos 2^_ a 2 b 2 + B 2 2 - S m 2 - 

2 V 2 



ttBF(B) 



o ° ■ 2 ° 2i2 , D 2 / n • 2 17 

2 — - sin - — a b + B 12 — sin - 



(CI) 

(C2) 
(C3) 

(C4) 
(C5) 



(C6) 



(C7) 



We have checked these formulae by performing the integrals with MAPLE and by direct numerical 
integration of Eqs. (B7) and (95) for a Gaussian potential. 
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